Fermionic propagators for 2D systems with singular interactions 
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We analyze the form of the fermionic propagator for 2D fermions interacting with massless over- 
damped bosons. Examples include a nematic and Ising ferromagnetic quantum-critical points, and 
fermions at a half-filled Landau level. Fermi liquid behavior in these systems is broken at criticality 
by a singular self-energy, but the Fermi surface remains well defined. These are strong-coupling 
problems with no expansion parameter other than the number of fermionic species, N. The two 
known limits, Ai" >> 1 and A'^ = show qualitatively different behavior of the fermionic propagator 
G{ek,(jj). In the first limit, G{ek,Lo) has a pole at some efc, in the other it is analytic. We analyze 
the crossover between the two limits. We show that the pole survives for all A*', but at small A'^ 
it only exists in a range 0{N^) near the mass shell. At larger distances from the mass shell, the 



system evolves and G{ek,Lu) becomes regular. At TV 
and G{ek,uj) becomes regular everywhere. 
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I. INTRODUCTION 

Physical properties of fermionic systems interacting 
with critical neutral fluctuations have been a focus 
of intense studies over the last several decades and 
yet remain a subject of advanced research. Exam- 
ples include fermions interacting with a gauge fieldfii^ 
a half-filled Landau level^^i^ and the system behavior at 
quantum-critical points (QCP) towards Ising-type ferro- 
magnetism'''^'^ (FM) and towards a nematic order ji^s^iiS 
The later case is an example of a Pomeranchuk-type 
Fermi surface instability of an isotropic Fermi liquictii. 
In all such systems, scattering of fermions by massless 
bosonic excitations leads to a non-analytic form of the 
fermionic self-energy S(/c,a;m)- Below the upper criti- 
cal dimension Dcr, 5](fcF,Wm) exceeds a bare term 
in the fermionic propagator, and the system develops a 
non-Fermi liquid behavior. At a critical point towards 
a nematic or an Ising ferromagnet, one-loop self-energy 

E(A:F,tfm) cc iLOm^ , and Dcr ~ 3 (Ref. [!!). The singu- 
lar behavior of the self-energy is, however, only in the 
frequency domain, the momentum dependence of 0) 
remains regular: S(A;,0) oc Cfc. As a consequence, the 
Fermi surface remains well defined at fcp as a locus of 
singular points of G(k, 0) despite that Landau quasipar- 
ticles do not exist. 



In 2D, one-loop self-energy becomes 



2/3 



It has 



long been the issue^i^i^ d^d^ whether form is the ex- 
act expression for a non-Fermi liquid fermionic propa- 
gator. The answer to this question is still lacking. On 
one hand, the two- loop and higher order self-energies also 



scale as 



2/3 



i.e., the exponent remains the same to all 



orders. On the other hand, higher-order terms are of the 
same order as one-loop self-energy, and it is a'priori un- 
clear what the sum of infinite series of terms yields. 

The way to treat such systems in a controlled way is 
to artificially extend them to N fermionic flavors and 
require that the interaction with a boson conserves the 
flavor. At large N , multi-loop self-energy terms ac- 



quire extra powers of (InA^/iV)^ and the series of uj^'^ 
terms converge. In this situation, self-energy is es- 
sentially determined by the one-loop term E(fc,a;„i) — 
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sgn(wm), where ujq is the internal energy 



scale. Accordingly, at 
below, 
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As a function of a complex _Fjv^oo(a^) has a simple 
pole at X = i. Because of the pole, real space/time prop- 
agator G{r, t) is long-ranged and decays by a power-law, 
as G{Y,t = 0) cx l/r^ and G{vFt » r) oc l/{t^). (see 
Appendix [A| . 

Another solvable limit is iV = 0. In this case the cur- 
vature of the 2D Fermi surface scales out, the system 
behavior becomes effectively one-dimensional and can be 
obtained by bosonization. The Green's function at A^ = 
has been obtained by loffe et al and Altshuler et al^ 
It is still given by Eq. (|1.1[) . but the functional form of 
F]s[=q{x) is fundamentally different: Fj^^q{x) is analytic 
in any finite region in the upper half-plane of x and be- 
comes singular only at x = oo. Because the pole is ab- 
sent, G'(r, t) is now short-ranged and at v^t ^ r scales as 
G{r, t) cx [l/t-yr] exp[— ar/t^/'^], where a is a dimensional 
prcfactor. 

Different behavior of Fn{x) at large N and at iV = 
raises the question which of the two forms (if any) de- 
scribes system behavior in the physical case of A^ — 1. 
Altshuler et al conjectured^ that the iV = case is spe- 
cial, and the system behavior at any A'^ ^ is qualita- 
tively the same as for large N , i.e., the pole in Fn persists 
for all TV and only vanishes at TV = 0. On the contrary, 
Fradkin and Lawler— argued that exponential behavior of 
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G(r, t) survives at finite N. One of us and Khveshchenko 
arguedi^ that the curvature of a 2D Fermi surface is rel- 
evant for any > 0, and G(r, t) decays by a power-law 
once > 0. However, the calculations in Ref. QJ^ are ap- 
proximate and did not yield the same G(r, t = 0) (x 
as at iV oo. 

In this paper, we analyze this issue in detail by per- 
forming loop expansion at small N. We find, in agree- 
ment with the conjecture by Altshuler et al— that the 
pole in the Green's function exists at any finite TV, and 
real-space, equal time Green's function G(r, t = 0) decays 
by a power-law, as The way how the pole disap- 

pears at TV ^ is, however, somewhat counter-intuitive. 
Naively, one could expect that the residue of the pole 
Zn gradually vanishes as ^ 0. We, however, found a 
different behavior: the residue Zn remains 0(1) at small 
A^, but the pole only exists in the range A < A^^, where 
A, introduced in (|1.2p below, is the distance from the 
pole. Outside this range, the Green's function is regular 
and the same as at A^ = 0. At A^ 0, the range collapses 
and the Green's function becomes regular even at A = 0. 

We present computational details below, but first sum- 
marize our rational. As our primary goal is to study what 
happens at small but finite A'^, we cannot use bosoniza- 
tion, which is only applicable at A^ = 0, and have to rely 
on the diagrammatic loop expansion. It is not guaran- 
teed a'priori that loop expansion is useful at small A^ as 
all terms in the series are generally of the same order, 
and it could be the case that the corrections are all reg- 
ular near the mass shell, but the prefactors are arranged 
such that infinite series of regular 0(1) corrections to the 
quasiparticle residue diverge at A^ = and destroy the 
pole. 

We, however, found that the actual situation is differ- 
ent, and the the pole disappears at A^ = and is still 
present at a finite N not because regular series diverge 
(at A' = 0) or almost diverge (at A^ > 0), but because 
of a peculiar singularity in the self-energy, whose form 
is different at A^ = and at finite A^. This singularity 
can be captured within the loop expansion. Specifically, 
we find that the expansion of the self-energy S near the 
mass shell, in powers of 



A = 1 



= 1 
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contains a universal term which is "non-perturbative" in 
the sense that it comes from fermions whose energies are 
of order ujA^^^, which are smaller than the external en- 
ergy u. This term appears at the two- loop order and at 
A^ = yields the non-analytic contribution to the self- 
energy which scales as A^/^. At a first glance, such term 
cannot eliminate a pole as it is smaller than Gq^ = A. 
However, this universal self-energy gets renormalized by 
logarithmically divergent vertex corrections, which, be- 
cause typical energies scale with A, are powers of In A. 
Series of such corrections exponentiate into the full ver- 
tex r cx A~° and modify the universal self-energy to 
A^/^r^ ~ The exponent a = 1 in the lead- 



ing logarithmical approximation (when only the highest 
power of the logarithm is kept at any order), but gets 
renormalized by 0(1) corrections beyond this approxi- 
mation. We cannot find a explicitly, but a comparison 
with bosonization implies that a = 3/4 in which case the 
fully renormalized self-energy tends to a constant value 
on the mass shell, i.e., the pole in G{k,uj) disappears. 

We next turn to A^ > 0. We found that universal 
contribution to the self-energy do exist in this case as 
well, but at the smallest A << A^^ it behaves as A^/N"^ 
(up to extra logarithms) Because typical internal ener- 
gies are still small, vertex corrections are again relevant, 
but now In A gets replaced by InA^^ for A << A^^, such 
that r cx l/N^^^. The full universal self-energy is then 
{A^/2/N^)r^ ~ (A/Af2)5/2^ ^j^g j.gg^l^^ ^jjg pQlg 

vives, at the smallest A, and its residue remains 0(1). 
However, {A/N^f^^ fo rm of the self-energy is only valid 
for A/N'^ . At larger deviations from mass shell, the self- 
energy approaches the same constant value as at A^ = 0. 
Generally, at small A'^, we have with logarithmic accu- 
racy. 
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Fn^o{x) = Fn^o (A) = 
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'A + g{A/N^)' 
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where g{y) cx y^/^ at y ^ 1, and g{y) ~ const at y ^ 1. 

To re-iterate, our key point is that there exist a uni- 
versal, non-analytic term in the loop expansion for the 
self-energy. At a finite A^, this term is cast into the scal- 
ing function of A/N'^, where A is the deviation from the 
mass shell. At the smallest A, it scales as A^/^ <c A, 
and the pole in G{k,ujm) then survives, with the residue 
Z + N = 0{\). AtA> A^^ the self-energy approaches a 
constant, and the full G(fc, LOm) looses the memory about 
the pole. At A^ = 0, the region where the pole exists 
vanishes, and the Green's function becomes regular even 
at A = 0. 

We also considered another example of a non-Fermi 
liquid behavior - the case of 2D fermions at the 
the half-filled {v = 1/2) Landau leveUJ^ii^ In this 
case, the fermionic self-energy is marginal at large TV, 
S(/c, ujm) — iAwml Incjml, whcrc A is a dimensionless cou- 
pling (see below). The fermionic propagator G(k, w,„) = 
~[l/ek)FN[-y), with y = Aa;m|lnw„|/efe then again 
has a pole at y = i. We solved the TV = limit 
by bosonization and found that F{){y) = e~*^, which 
is obviously a regular function along imaginary y axis. 
We then performed the small TV analysis and found 
the behavior which is similar but not equivalent to the 
previous case. Namely, at the distance from the pole 
A = 1 + iefc/(AtJ„i| Ino^ml) ^ TV^, the pole still exists at 
any finite TV, but its residue now scales as Zm oc TV^. At 
larger deviations from the mass shell, the universal self- 
energy approaches a constant, and the system looses the 
memory about the pole. At TV = 0, the range where the 



pole exists collapses, and the Green's function becomes 
regular for all A. 

The paper is organized as follows. In Section II we 
introduce the model. In Section III we consider a 2D 
system at a Pomeranchuk QCP towards a nematic order. 
We present the results for large N and briefly review 
bosonization results for = 0. We then discuss universal 
terms in the loop expansion at A'^ = and at a finite N. 
In Section IV we present the same consideration for 2D 
electrons at the half-filled Landau level. In Section V 
we present the Conclusions. Some technical details are 
presented in the Appendices. 



II. THE MODEL 

We consider 2D fermions with a circular Fermi sur- 
face and dispersion e^. We assume that fermions inter- 
act at low energies by exchanging collective excitations 
with the static propagator x{q) — Xo/q^'^^- For nematic 
and Ising ferromagnetic QCP x = Ij^ili^i^ili^iiS for the 
half-filled Landau level with unscreened Coulomb inter- 
action, a; = (Refs. l2 Hl"5lfT6l ) . We will only consider 
interactions with charge fiuctuations. Interactions with 
gapless spin fiuctuations require separate treatment .A§ 
The static xil) is predominantly created by high-energy 
fermions and is an input for the low-energy model The 
Hamiltonian of the model is given by 

k,a q 

+ 9 ^4,a^"/3Cfe+«,/3 (2.1) 
k,q 

where the first term is the kinetic energy of fermions, 
the second term is the potential energy of collective ex- 
citations described by and the third term describes 
the interaction between fermions and collective modes. 
The coupling g generally depends on momentum but 
can be approximated by a constant at small momenta. 
Pai3 — <^af3 Ising ferromagnet and Sap for a nematic 

transition. The results are equivalent in both cases, and 
we will only consider a nematic transition. 

Near a particular kp point, which for definiteness we 
direct along x, 

ek = [(fcx - fcp) -I- kl/{2k^)] . (2.2) 

The second term is due to the curvature of the 2D Fermi 
surface. 

The interaction g appears in the perturbation theory 
only in even powers, in a combination g — g^Xo- The 
dimension of g is inverse mass 1/m, the dimension of xo 
is mk^, hence the dimension of g is energy. The model 
has a natural dimensionless parameter A = g/E-p, where 
Ep is the Fermi energy which we assume to be of the 
same order as the fermionic bandwidth. We assume, like 
in earlier worksj^i^ that A << 1. This condition implies 
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FIG. 1: The one- loop self-energy diagram. 

that interaction does not take the system out of the low- 
energy domain, i.e., low-energy behavior is well separated 
from the system behavior at energies compared to the 
bandwidth. 

As it is customary for the problems in which fermions 
interact with their own collective modes, collective exci- 
tations become Landau overdamped due to interaction 
with fermions, and the full dynamic susceptibility of the 
<& field becomes 



where 7 = gkp/{Trvl) (see Refs.HH). The Landau damp- 
ing can be included into the theory already at the bare 
level, all one has to do is to change from a Hamilto- 
nian description to a description in terms of an effective 
actiouiii 

The model can be extended to A^ ^ 1 fermionic fia- 
vors by adding a fiavor index to fermions and keeping 
flavor index intact in the interaction with the 4> field. 
This extension allows one to consider large A'^ and small 
A^. At large A^, fermionic damping 7 scales as N and 
is large. Collective excitations then become slow modes, 
and their effect on fermions becomes small in 1/N by 
Migdal theorem. This is the limit where a direct pertur- 
bative treatment is applicable. In the opposite limit of 
small A^, the momenta which mostly contribute to the 
fermionic self-energy are of order A'^, such that the ky 
term in Eq. (|2.2p is small by A^, and the curvature of the 
2D Fermi surface becomes a small perturbation;^i2, With- 
out the ky term, the fermionic dispersion becomes purely 
one-dimensional, and the self-energy can be found by a 
bosonization technique. As we said, our key goal is to an- 
alyze the crossover between solvable large A^ and A^ = 
limits with the aim to understand system behavior for 
the physical case of A^ = 1 . 

Below we consider separately the case of short-range 
(screened) interaction for which xil) l/?^ (this is the 
case of fermions interacting with quantum-critical collec- 
tive excitations), and the case of unscreened long-range 
interaction (electrons at the half-filled Landau level) for 
which x{q) oc l/\q\. 

III. SHORT RANGE INTERACTION, A 
QUANTUM-CRITICAL POINT 

The limiting cases N ^ 00 and A^ = have been 
studied before. We briefly review the existing results and 



go one step further in the 1/7V expansion for large N. 
We then present and discuss our results for the universal 
terms in the loop expansion near the mass shell. 



A. Fermionic propagator at A'^ 3> 1 

At large N, the results for the fermionic propagator 
can be obtained by expanding in the number of loops. 
Each extra order brings extra smallness in 1/JV. Explicit 
calculations show that the parameter for the loop ex- 
pansion is actually In'^ N/N'^ , which is even smaller than 

The one-loop self-energy diagram is shown in Fig. [TJ 
The k dependence of the one-loop self-energy is regular, 
with a prefactor 0(A) << 1, and we neglect it. The 



frequency dependence of Ei(fc,ct;,„) is uj: 



2/3. 
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The prefactor ujq formally contains 1/iV, but it can be 
absorbed into the renormalization of the Fermi velocity. 
All higher-order diagrams contain the same combination 



2/3 1/3 



and we will just consider loq as a normaization 



factor for frequency. 

The a;2/^ dependence of E implies that quasiparticles 
are not sharply defined (along real frequency axis, ReT, 
and /mS arc of the same order uj'^^'^), i.e., the system be- 
havior is not a Fermi liquid. At the same time, S(fc, ujm) 
still vanishes at zero frequency, i.e., the Fermi surface re- 
mains well defined. Substituting the one-loop self-energy 
into G(k, and neglecting LUm in comparison with Si 
we reproduce (jl.ip . 

The relevant two-loop self-energy diagram is shown 
in Fig. [2^1.. Earlier estimates show that, when = 
0, S2(fcp,tJ™) oc iojK^ujl^\\nN/Ny (Refs. |2|l5|). We 
found that S2(fc,a'm) is also a non-trivial function of 
ujm^ujg^^ /ek- For our purposes, we need the expres- 
sion for 'S2{uJrmk) near the mass shell A = 0, where 

A = 1 + iek/itUm^ul^^). We found (see Appendix B) 
that, to order A^, this expansion is regular and 

2 



,3 1/3/lniVy 



(3.2) 



0.52-F — A-h0(A2) 
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Adding E2 to Ei, substituting into the Green's func- 
tion and casting G into the form of Eq. , we find 

Zn 



Fn{x) 



Dn + ix 



^" -Fine (-^) ; 



(3.3) 



where Fine is a regular function near the mass shell x — 
—iDf^, and 



/lniV\' 



Dn = 1 -0.52 
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FIG. 2: (a) The relevant two-loop self-energy diagrams. The 
universal contributions come from oji,ij — UJ2 ~ ujA'^^^ <^ uj 
and 102,1^ — i^i ~ 

^3/2 

LO <^ u}. (b) Relevant vertex corrections 
to the two loop diagram whose universal contribution comes 
from LUi,uj — UJ2 ~ uiA'^^^. For the second universal contribu- 
tion, from uj2,i^ — i^i ~ A'^''^, the corrections to the other two 
vertices are relevant, (c) The full universal self-energy with 
fully renormalized vertices. 



We see that, at this level of consideration, correc- 
tions lead to three effects: the pole in G{ek,x) acquires 
the residue < 1, the location of the pole along the 
imaginary x axis shifts to a somewhat different x = iD^ 
and the Green's function acquires an incoherent part 
which is regular near the pole. Perturbation theory in 
1/iV is perfectly well defined, and the pole in Fn{x) is 
surely present at large N. 



B. Fermionic propagator at A'^ = 

In the limit = 0, the curvature of the Fermi surface 
scales out and the original 2D problem maps onto an ef- 
fective ID theory of electrons with retarded interaction 
(Refs. [2113) ■ This allows one to compute fermion prop- 
agator by using ID bosonization technique. The result 
for G(k, LUm) can again be cast into the form of Eq. p.ip . 
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FIG. 3: The function FM=o{iy)- Dashed line represents small 
y asymptotics (see text). 



but now^ii^ 

F^^oix) = ^exp{(-l)5/4r3/2(5/3)x3/2} 

exp{i(r(5/3)a;y)3/2} 



2 ' 



dy 



An ./o + iy-l 



(3.5) 



The analysis of Fn^q{x) is presented in Appendix A. The 
key point is that Fn^q{x) is a regular function in any 
finite region in the upper half-plane, and, in particular, 
along imaginary x — iy,y > 0, where there was a pole at 
large N. At small x, Fn=(){x) expands in regular powers 
of x: 



FN=oix = iy) = l + y + 



^^rHWS)y' + 0{y% (3.6) 



At large x, Fj^=o(x = iy) grows exponentially with y, 
but still remains finite for all finite y. We plot FN=o{'iy) 



in Fig. [3 



C. Nonanalytic terms in the self-energy 

Our key interest is how the pole disappears between 
= and >> 1. One possibility might have been 
that the pole drifts to higher x with decreasing A^, and 
disappears at a; = ioo at N — 0. That would be consis- 
tent with the divergence of F]\[=o{iy) at y = oo. However, 
the expansion in 1 /N at large A^ shows that the pole actu- 
ally drifts to smaller x = iy with decreasing A^. Another 
possibility would be a behavior in which the residue Zn 
of the pole gradually disappears at A — > 0, i.e., regular 
corrections to Zn (the ones which are of order In^ N/N"^ 
at large A^) make Zq = 0. We, however, didn't find a 
self-consistent solution for Zjv and D^r at small A^ (this 
search required some computational efforts). 

Below we assume that regular corrections to fermionic 
propagator leave Zm and Dm finite at A^ — > 0, and the 
pole is destroyed by "universal" non-perturbative terms 
in the loop expansion of the self-energy near the mass 
shell. These terms appear at every order in the loop ex- 
pansion, are non-analytic, and come from fermions with 
the lowest energies. We first identify these terms at 
A^ = 0, and then show how they get modified at finite A^. 



1. N=0 

The two-loop diagram is presented in Fig. [5] contains 
momentum and frequency integrals. At A^ = 0, the cur- 
vature term that couples the two disappears (see Eq. lB2p . 
The momentum and frequency integration in the self- 
energy then factorize, and the momentum integrals are 
straightforwardly evaluated leaving only non-trivial fre- 
quency integrals. The two- loop self-energy becomes 



E.(fc,.„) . l^-j (^^j ^L,^ 

1 

A + [(1 - r!i)2/3 + (1 _ r!2)2/3 + (r!i + ^2 - 1)2/3 _ f] ' 



where ili and ^2 are two internal frequencies normalized 
by ujrm and, we remind, A = 1 -I- ick/ [iio'U.'^ujy^). Note 
that the coupling g is fully absorbed into wq, i-e., for a 
generic A = 0(1), T,2{k,u}m) (as well as self-energies of 
higher-loop order) is of order , with a prefactor 

of order one. 



Expanding Eq. p.7p to first order in A, we only obtain 
regular contributions to Zq and to Dq. These contribu- 
tions are perturbative in the sense that typical internal 
are of order one, much larger than A <C 1 over 
which one expands. However, this doesn't go beyond 
first order - a formal expansion to order A^ yields a di- 
vergent prefactor. On more careful look, we found that 
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the next term beyond A is actually A^/^. This term 
comes from the regions ili << 1, 1 — — ^2 << 1 
and << 1, 1 — ill — ^1 << 1- The contributions 
from both regions are equal, and we only focus on the 
first one, for which the denominator in Eq. (j3.7p becomes 



A + (11^)2/3 + - fl^)2/3. Rescaling fli = A^^^x, 
n'2 — A^/^y and multiplying the self-energy by a factor 
of 2 to account for the two regions, we re-write Eq. p.7p 

as 



, A3/2 

\p ^ o- 2/3 1/3 ^ 



(3.8) 



The largest contribution to this integral is confined to 
the upper limit of the integration over x. This is a reg- 
ular, perturbative term (internal frequency (rii)^/^ — 
Aa;2/3 /\y However, we also found a contribution 
which is confined to x — 0(1), i.e., to internal energies 
of order A^/^. For the integral in (|3.8p this contribution 
can be easily singled out - it is given by 



+ 



(3.9) 



^ o. 2/3 i/3A=^/2 dx 



Q / 7-1/3 , 
Jo ^ ' Jo 

1 



(1 + y2/3 + (a; _ 2^)2/3)(y2/3 _^ (2; _ y)2/3)2 ' 

Evaluating the integral numerically, we obtain 

S2(A:, uJm) = 2.69 zc^^/^w^'A^/^. (3.10) 

At the first glance, this term is smaller than 0(A) and 
is irrelevant to the issue of pole disappearance. Let's, 
however, continue our analysis and verify whether there 
are singular renormalizations of £2(^1 ^m) from three- 
loop and higher-order diagrams. We found that such 
renormalizations do exist and come from vertex correc- 
tions. In particular, for the universal contribution to 
S2 from ill << 1, 1 — ri2 = ^'2 << 1, singular renor- 
malization comes from corrections to the two spin-boson 
vertices which involve bosonic propagator with the fre- 
quency Q2 ■ Each vertex correction contains a block made 
out of two extra Green's functions and one interaction 
line. We verified that the corrections to the two vertices 
involving the propagator with Q2 are the same, and each 
yields, after integrating over internal 
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(27r)2«, y^^^ Si(i7) J r 

2 dz _2 1 



z 



(3.11) 



We see that vertex correction is logarithmic and large, 
if ill ^ 1- Typical fli for the universal term in E2 
are of order A^/^, hence the leading vertex correction is 
ln(l/A). Observe that the prefactor is just a number - 



the coupling g is canceled out by {'-fujoy^^ oc g. Evaluat- 
ing higher-order corrections in the leading logarithmical 
approximation (in which we keep only the highest power 
of the logarithm at any order) , we find that the full vertex 
r becomes (see Fig. [2j;) 



r 
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(3.12) 



This r is the solution of the differential RG equation 

1 



dr 

dL 



r, L = In ■ 



(3.13) 



The implication is that the vertex correction can be 
treated within RG technique, i.e., that the system is 
renormalizable . 

Eq. (|3.12p is not exact, however, as there exist extra 
terms of order ln(l/A) from second-order and higher- 
order vertex corrections. Such terms come with prefac- 
tors 0(1), and modify the exponent in (I3.12p to 



F cx 



(3.14) 



where a = 0(1). Because the exponent a has contri- 
butions from all orders, we cannot compute it explicitly 
within the loop expansion, nor verify explicitly that log- 
arithmic series still exponentiate beyond the leading log- 
arithmic approximation (i.e., that RG procedure leading 
to ()3.14p is still valid). At the same time, the situa- 
tion here is not different from a variety of other prob- 
lems where RG treatment has been adopted based on 
the leading logarithmic series ^iiiii^ and we assume that 
Eq. (|3.14p is valid without further reasoning. 

Combining Eq. ([SH]) with Eq. ([STO]) . we find that the 
universal part of the self-energy, dressed by logarithmic 
vertex renormalizations, behaves near the mass shell as 



(3.15) 



To agree with the bosonization formula, we have to set 
a = 3/4, then Eun™ becomes A- independent, and the 
full Green's function behaves near the former mass shell 
as G~^{k,ujm) — iuhn'ujy^ (A -|- BA + ...) where A and 
B are constants. 

We now re-evaluate universal terms at 7^ and see 
how the universal self-energy gets modified at a finite A^. 



2. Finite N 

At a finite N, momentum and frequency integrals do 
not decouple, and the computations become more in- 
volved. Still, the two momentum integrals can be evalu- 
ated exactly, at any A^, so the remaining task is to prop- 
erly estimate the frequency integrals. As at A^ = 0, we 
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expand in A and verify whether or not this expansion 
contains the non-perturbative, universal term. We find 
that such term is again present, but at a finite N and 
A < N^, scales as A^/^ rather than A'^/^. This term 
comes from the same range of internal frequencies as at 
A'^ = - one of the internal frequencies is small, another 
is close to external lo. Using the same notations as at 
N = 0, we obtain for the contribution from such region, 
with logarithmic accuracy 

^5/2 



^ 

■ N 
1 



2/3,„l/3. 



In 



47r2iV2 

/ ( 1 + 2/2/3 + _ 2^)2/3 
2.1/3 



(3.16) 



■In 



I + y2/3 + (2; _ ?;)2/3^ 

The universal term (the one which does not depend on 
the upper limit of integration) comes from the cross- 
product of the two logarithms, and the universal part 
of the two-loop self-energy is 



E2(fc,c.„0=^c.2/3^y^A5/2 
The prefactor J is given by 

J = — / dx 

where 

V'(a;, z) = \vl\\ 



J 



87r27V2 



In 



Ar2 
A" 



dziplx, z) 



1 



In 1 



^2/3/. 
1 



1 



X 



a;2/3/, 
1 



X' 



.2/3j 
1 

2 



(3.17) 



(3.18) 



(3.19) 



and fz = z^^"^ + (1 — 2;)^/3^ Evaluating the integral, we 
obtain J « 0.96. 

The next step is to include vertex renormalizations 
from higher-order diagrams, i.e., renormalize the univer- 
sal part of E2 in p.l7p into "ili^univ = S2 * T^- We found 
that multi-loop contributions to F are again logarithmic, 
but now the lower limit of the logarithm is set by the 
largest of A and N"^ . The fuU result is then V cx l/A^/^ 
for A > 7V2 and V cx 1/7V3/2 for A << Ar2 Combining 
this with (|3.17p . we obtain 



2/3 1/3 { ^\ 



(3.20) 



where g(ij « 1) cx y^/2 and g{y >> 1) — const. We 
neglected In A/iV2 term because F is only known to the 
accuracy that neglects possible additional logarithmical 
factors. 

Substituting the full T^univ into the Green's function, 
we reproduce Eq. ([TT 



G{k,UJm) = 

Fn^o (A) 



--Fjv^o(A); 

1 

-A + g(A/A^2)- 



(3.21) 



We see that the pole survives, at the smallest A for any 
nonzero N and, moreover, its residue remains Zjq = 
0(1). However, at distances from the pole larger than 
N"^ , the system evolves, and its behavior becomes essen- 
tially the same as at TV = 0. 



IV. HALF-FILLED LANDAU LEVEL 

In this section we study different example of a non- 
Fermi liquid behavior - 2D fermions at half-filled Landau 
level with unscreen Coulomb interactioni^ i^^'^^ The ef- 
fective low-energy theory is described by Eq. (|2.ip with 
massless bosonic propagator 



muq + 7 \uj/q\ 



(4.1) 



Here u is the effective velocity [u = /{Sire) with e the 
dielectric constant of a host semiconductor], and 7 = 
gk-p/{TTv1) is the same as in the Li;2/3 problem. 

We start by presenting the results for the limiting cases 
N ^ 00 and = 0. Large N limit is studied in 1/A^ 
expansion while iV = is obtained via ID bosonization. 
We then discuss our results for the universal terms in the 
loop expansion near the mass-shell. 



A. Fermionic propagator at >> 1 

The one-loop self-energy has been calculated in Ref. [3. 
It has the form 



Si(w„) = iX ujm In 



Jo_ 

\uJrn 



(4.2) 



where X = g/ (^Air^ukp) is a dimensionless coupling and 
eo is an energy scale defined by eq — (^tt/N) [E^/g) uk^. 
The momentum dependent part of S is regular and we 
neglect it. The ujlogui form of the self-energy Eq. (|4.2p 
implies that the system exhibits a marginal non-Fermi 
liquid behavior, however the Fermi surface remains well 
defined. 

As in the previous case, higher order self-energy dia- 
grams are parametrically small in (lnA^/A^)2. Near the 
mass-shell A = 0, where 



A = 1 



Xcurn ln(eo/a;m)' 



(4.3) 



we found at two-loop order (see Appendix [C]) that 
Y,2ik,i-u„i) is linear in A 

^2ik,ujm) = iXujm'^'a{eo/uj,„) ^ (^-^^] . (4.4) 



24 V iV 



We see that, as for the ^2/3 problem, multi-loop self- 
energy terms acquire extra powers of (lnA^/iV)2. Sub- 
stituting one- and two-loop self energies into the Green's 
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function and rewriting it as in Eqs. (II. ip . (|3.3p . we obtain 
G(k,w,„) = - (l/ek) Fn^-v), where y = (Si/iefe), and 
pN^ooiy) is again given by Eq. (|3.3p . but now 



1 

24 



In TV 



D 



N 



(4.5) 



We see that 1/N corrections for the Landau level case 
are qualitatively similar to the case of a nematic QCP - 
in both cases, the pole in G(k, cjm) viewed as a function 
of e/j. survives, but acquires a residue Zj^ < 1. 
In real space, we have (see Appendix IC II) 



G(r) = Go(r) 



1 



, vpt < rlnr/ro, 



_ln(r/ro)_ 

G(r,i) «Go(r,i), t-^t > r ln(r/ro), (4.6) 



where Go(r, i) cx (r + Wpt) ^^kp/r is a propagator of a 
free fermion, and tq = fp/eo- 



B. Low-energy effective theory for = 0: 
Fermionic propagator 

As in the previous case, at = the curvature of 
the Fermi surface becomes unimportant, the motion of 
fermions becomes essentially one-dimensional, and the 
fermionic propagator can be obtained by mapping the 
original 2D problem to an effective ID theoryi^ii^ For a 
half-filled Landau level the effective action has the form 



S ^ J dndk (4.7) 



where p^^k = '0° k density operator, and the 

replica index a has s values. As usual, we take the limit 
s in order to avoid generating fermionic loops in 
perturbation theory. 

We solve the model Eq. (|4.7p by ID bosonization. Fol- 
lowing the standard steps we find that 



G^^(r,T) 



27r(r — ivj^r) 



exp 



-A 



r In |r — 



(r — iUpT) 



(4.8) 



where t is Matsubara time. A Fourier transform of 
Eq. (|4.8p yields a simple scaling form for G(k, cj„i): 



G(k,Lj„i) = — -Fn=o 
Fn=o{x) = e-". 



-Xu/m In {eo/\uim\) 



(4.9) 



We verified that the loop expansion reproduces the small 
X expansion of Fi^^q{x). Obviously, Fi^^q{x) has no 
poles along x — iy where ^V=o(*y) — e^. We verified 
that the singularity at y — * cxd is responsible for the ex- 
ponential behavior of G{r,t) in ([47 



C. Nonanalytic terms in the self-energy 

1. N = 

The relevant two-loop self-energy is again given by the 
diagram in Fig. [5^. We compute it in Appendix [C] for 
arbitrary N. As in the previous case, consider first N — 
0. We have 



r-l rl 
1^2{k,UJm) = iX UJ,n I dQi / dfl2 







1-Sli 



In [eo/(w„irii)] ln[eo/{uj„,,Q2)] 



A In (eo/wm) + {(1 - r!i) In [^] + (1 - ^2) In [^] + (fti + ^2-1) In 



(4.10) 



where fli and ^^2 are the two internal energies normalized 
by ujm- In the previous case, S2 was regular to first order 
in A. This time, expanding Eq. (|4.10p to first order in 
A we found that the prefactor diverges logarithmically. 
This suggests that the term next to the constant is al- 
ready nonanalytic. We explicitly verified that the nonan- 
alyticity emerges in Eq. (|4.10p from the regions fli 1, 
572 ~ 1 and ^l2 ^ 1, i^i ^ 1- Within the logarithmic ac- 
curacy, the nonanalytic contribution to E2 can then be 



cast into the form 



Y.2(k,uj„i) = -2iXujm A In (eo/u;„i) 
"1 dfli 

ln(.o/c>„„) Jlllnl/fil 
|1„A| 

-2iX u;,nA In^ (eo/w„0 Inln (^j . (4.11) 
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As before, higher-order diagrams modify Eq. (|4.1ip by 
adding vertex corrections. A building block for a vertex 
correction is again the product of two fermionic Green's 
functions and one interaction line. Evaluating this block 
we find that vertex corrections are again logarithmical. 
Using further the fact that typical dimensionless external 
energy for the correction to the vertex involving a boson 
with Q2 is of order Oi , we obtain the renormalized vertex 
r in the form 



f = 1 + ln(l/f7i) 



(4.12) 



Evaluating further higher-order vertex corrections, we 
obtain in the leading logarithmical approximation 



(4.13) 



As in the previous case, the exponent is modified by 
higher-order corrections to 



(4.14) 



A simple experimentation shows that the full universal 
self-energy T,univ = ^2^^ agrees with bosonization if we 
set b = 1/2. Indeed, in this case, we have, combining 
Eqs. and ([iTT|) 



umv 
1 



-ic\ Lu„i^\v? (ea/LOm) x 
dVLi 



(4.15) 
icXujjn ln(eo/t^m) , 



where c — 0(1). We see that the fully renormalized 
tends to a finite value on the mass shell, i.e., 



G ^(A) = iAti^m In (eo/tiim) A -f c 



(4.16) 



This behavior is consistent with the bosonization for- 
mula, Eq. (|4.9p which yields oc (1 -I- A) near the 
mass shell. The computation of the constant c is beyond 
the scope of our analysis. 



2. Finite N 

We now check how the expression for the self-energy 
is modified at finite N . Just like at = 0, we expand 
the two-loop self-energy in powers of A and extract the 
universal term in the prefactor. We find that, to logarith- 
mical accuracy, the linear in A term is already universal, 
like at A^ = 0, and the only difference between A^ = and 
N > cases is in lower limit of the frequency integration 
over fJi: at a finite A^, instead of A, as in (|4.1ip . we now 
have A^^. Accordingly, instead of ( (|4.1ip . we now have 



'iXujrn lii^ (co/^m) Alnln(l/A^). 



As the next step, we include vertex corrections. They 
are still logarithmical and exponentiate, as in (I4.14p . but 
now typical fli arc of order A^^, hence 

2 A 

cx -lA cj,„ In (eo/t^m) (4.18) 

We neglected In | In A^| factor as the renormalized vertex 
is only known to this accuracy. We see that Tiuniv is still 
linear in A, i.e., the pole in the Green's function survives 
in the range of order A^^. This is similar to w^^'^ case. 
Unlike that case, however, the residue of the pole scales 
as A^^ and vanishes when N = 0. Substituting 'Suniv into 
the Green's function and rewriting the result in terms of 
F/v, as in H.9L we obtain near the mass shell 



G(k, w„ 



— Fn- 



.(A) 



Fn^o A] =F 



A 

7V2 



(4.19) 



(4.17) 



where F{y ^ 1) cx y, F{y 3> 1) = const. 

The difference in Zjv between the QCP and the Lan- 
dau level cases is due to the fact that in the case of a 
Landau level the universal self-energy already emerges in 
the loop expansion already at order A while at a nematic 
QCP it appears at order A'^/^. There may also be ad- 
ditional differences between the scaling functions -Fjv-»o 
and Fjv^o due to extra logarithmic or doubly logarith- 
mic factors, but, as we said, these factors are beyond the 
slope of our paper. 



V. CONCLUSIONS 

In this paper we considered two examples of 2D 
fermions coupled to critical overdamped bosonic fields. 
One describes a nematic and an Ising ferromagnetic 
quantum critical points, other describes a half-filled Lan- 
dau level. For both cases, the low-energy physics is 
governed by the interaction between fermions and col- 
lective neutral excitations with the static propagator 
x{q) = Xo/q^~^^- For a nematic and ferromagnetic QCP 
cc = 1, for a half-filled Landau level with unscreened 
Coulomb interaction, x = 0. In both cases, quantum fluc- 
tuations destroy a coherent Fermi liquid behavior down 
to the lowest energies, but leave the Fermi surface in- 
tact. The issue we addressed is what is the form of the 
fermionic propagaror, and whether it has a pole as a func- 
tion of the dispersion e^. 

The low energy properties of these systems can be 
studied in a controllable way by extending the theory 
to A^ >> 1 fermionic flavors and assuming that the 
interaction with neutral excitations conserves the fla- 
vor. At large A^, self-energy is perturbative in 1/A^, and 
the pole in the Green's function survives despite that 
^(fcF,^) cx w^/'^ for X = 1 and 'S{kp,uj) oc wlogw for 
X = 0. 1/N corrections only affect the residue of the 
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pole and make Ziq < 1. Existence of the pole implies 
that the propagator in real space/time is long ranged and 
decays by a power-law, e.g., as G{r ^ v-pt) cx l/r^ and 
G{r < Vpt) cx l/ty^ for x = 1. The other limit N = 0, 
has been solved using ID bosonization, and the result 
is that for both x ^ 1 and x — 0, fermionic propaga- 
tor does not have a pole, and the propagator G'(r, t) is 
short-ranged. 

The issue we addressed is at what N the pole disap- 
pears. This is essential as the physical case = 1 is "in 
between" the two limits. We performed the loop expan- 
sion for the self-energy both at TV = and at a finite 
but small N. At N — 0, we identified universal, non- 
analytic contributions to S which destroy the pole and 
make fermionic propagator regular near the former mass 
shell. At small but nonzero N, we found that singular, 
universal terms in the self-energy still exist, but they do 
not destroy the pole in the range of order 7V^ around 
the mass shell. At larger deviations from the mass shell, 
fermionic propagator recovers the same regular form as 
at = 0. For a nematic QCP, the residue of the pole 
Zjv remains 0(1) in this range, while for the case of a 
half-filled Landau level, the residue of the pole scales as 
Zn cx N^. 

The result that the pole in G{k,uJrn) exists for any 
> and only vanishes at = agrees with the con- 
jecture by Altshuler et The key result of our work is 
the identification of peculiar, universal terms in the loop 
expansion of the self-energy, which are responsible for the 
destruction of the pole at A^ = 0. At A^ > 0, these terms 
get modified - they do not destroy the pole but reduce 
the width where the pole exists. 
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APPENDIX A: ANALYSIS OF THE FERMION 
PROPAGATOR AT A QCP AT Af ^ oo AND N ^0 

In this Appendix we analyze the behavior of the 
fermion propagator 

at a nematic QCP at A^ >> 1 and TV = 0. The goal 
of this analysis is to demonstrate that the form of long- 
distance behavior of G{r, r) , where r is the Matsubara 
time, is qualitatively different depending on whether or 
not G{\i,ujm) has a pole as a function of e^. 



We use the Matsubara frequency form of G{ek,ujm) = 

/ 2 /3 1 /3 \ 

— {l/€k)FN [~^m /cfc)- Introducing new vari- 
ables, X = -(wm^tjQ^^)/efc, y = {ekr/vF) and a = 
I^WF^^r/t^y^r'^/^^ we re-write Eq. (|Aip as 

X cos(fc,r + y-J)F^(.)e-(-i)^^^-^^V/^". 

(A2) 

1. Limit N oo 

In this limit Fn^oo^x) = 1/(1-1- ix). Substituting this 
into Eq. (|A2p . we obtain 

■/ -I \3/2 3/2 3/2 

g-H-l) ^ X ^ y ^ a , 

" X^x COS (fc.r + ,--]. (A3) 

We now show that the pole sX x — —i determines long- 
distance behavior of G'(r, i). For definiteness, we focus on 
the case a « 1, i.e., Wf^^t << o^Q^^r"^/^. In this limit, 
the exponential factor in (|A3[) can be set to unity and the 
integrals decouple. Integration over y is straightforward 
and yields 




We added the factor e '^'1^1 for regularization. The sub- 
sequent integration over x is also straightforward. 

/■oo I 11/2 roo 1/2 

/ dx^ = 2 / dx— — ^ = \/2^. (A5) 

J-oo 1 + «a; Jo 1 + a;-^ 

We also evaluated the integral in Eq. (jA5[) by extending 
the integration into the complex plane of x and found 
that the result comes from the pole at a; = i, the branch 
cut contribution in (IA5[) is cancelled out. Substituting 
Eqs. ((X4| and (IX5|l into Eq. (|X3)) we find that the fermion 
propagator G(r, r) cx l/r^ at vpT <C r^^'^lrl/'^, where 
To = vp/ujQ. We emphasize again that the power-law 
decay is the consequence of the pole in G(k, ojm)- 

For completeness, we also present the result for G(r, r) 
for a ^ 1 . In this case we obtained 

G(r,r)^lf^),/^/C(a), (A6) 

where 

IC{a)\ ^^cosffcpr-^V (A7) 

^ '\a^oo 3q, V 4/ ^ ' 
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Using a cx r/r^/"^ we find that in this hmit the the integration contour into the complex plane. As for 
real space/time Green's function behaves as G(r, r) c>c large iV, the contribution from the branch cut is canceled 



1/ {-^v-pt). Combining the two results, we have 



G(r,r) 



CX 



, (A8) 



out and in the absence of the pole G(r, r) only comes from 
the integral over a semi-circle with infinite radius. Using 
parametrization 



3/2 _ p3/2 



R-^' e^p {iip} , TT < < 27r, 



(A14) 



or, in terms of the frcc-fcrmion propagator Go(r, t) 



Go(r, r)yV^, VpT < r^/Vr^^ 



with R ^ oo, substituting into (|A2ll the large- a; asymp- 
tote of Fn=o 



(A9) 



Fn=o{x) = - exp 



jTv/4+i^-p3/2 



r^/"(5/3)i?3/2 , (A15) 



and introducing yi — a^^^y, we obtain 



2. Limit iV = 



1 



The expression for G{k,u!m) in this case is given in 
Ref. 0, see Eq. The function Fn =oix) in (|3.5p can 

be represented as 

FN=a{x) = cosh\e'''/'^r^/'^{5/3)x^^^]-ixTi{x) 

- x2^r3(5/3).F2(x), (AlO) 

where J-'i and are expressed in terms of hypergeomet- 
ric functions 

Ti{x) = 1F2 [1; 5/6, 4/3; (l/4)ir3(5/3) x^] , 

T2{x) = 1F2 [1; 7/6, 5/3; {l/A)iT^{5/3) x^] . (All) 

At small arguments, J^i{x) and J-2{x) are expandable in 
:Fi{x) = 1 + -^ [r{5/3)x]' ^ [r(5/3)x]^ + O [^x'] , 



G(r,r) = . 



X cos I k, 



27r3/2 
yi/aV3 



X - exp 
2 



(A16) 



Integrating over the angular variable we obtain 
i?3/2^"'d{e^^}exp|i?3/2e^^ 

"_,y3/2 _ g^V4r3/2(5/3)- 
",2^3/2 ^g„/4p3/2 (5/3)1 ^ 



= 27r(5 



(A17) 



^.(x) = l + A.[r(5/3).]^-^ 



[r(5/3)x]*' 4- O 



where 5{...) is the (5-function. The integration over yi is 
then straightforward and yields 



By contrast, at large |a;| ^ 1, Fn=o{x) diverges as 
3 



cos I 



Fm=o{x) 



2 exp 



g<V4r3/2 (5/3)^3/2 

3i^/3 2 



ar^ \7r^/2 

2/3 



yi/a 



(A13) 



Attx 9r3(5/3)a;2 



O [1/x^] 



: exp 



r(5/3) 



-/ dy,\y,\'/'x 
)5[z2/?/' + e^^/4p3/2(5/3) 

(A18) 



(i)2/3Q,2/3 

Re-expressing a in terms of r and we finally obtain 



r(5/3)r 



1/3 r- 1 



2/3 



(A19) 



Still, along x — iy F/v=o (iy) as a regular function of ?/ for 

all finite y, i.e., there is no pole. The plot of FM=o[iy) G{y^t) oc — —\j — exp| 

versus y is presented in Fig. [3l 

Below we show that the absence of the pole eliminates 

a power-law decay of G(r, r), while the divergence of We see that G(r, t) is now exponential. The full 

Fn^q(x) at infinity gives rise to the exponential decay bosonization result contains (wpT — r) instead of WpT 



of the fermionic propagator. 



in denominator of (|A19|) (see Ref. [l3) . To reproduce it. 



The easiest way to evaluate the integral over x in we would need more complex form of G(A:, Wm) than the 
Eq. (|A2p with Fn=o{x) given by Eq. (|A10[) is to deform one we borrowed from Ref. [3- 
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APPENDIX B: TWO-LOOP CONTRIBUTION TO THE SELF-ENERGY AT A QCP. 

In this Appendix, we show the details of the calculations of the two- loop and three-loop self-energy diagrams. 



1. Two-loop diagram 

The diagram is presented in Fig. [2l In analytic form 
4 f duji dqi f duj2 dq2 



x(qi, t^i)x(q2, t^2)G'(k + qi, - cji) 



(27r)3 J (27r)3 
X G(k-f q2,cj„i -W2)G(k + qi +q2,t^„i - c^i -6^2) 



(Bl) 



9 



(27r)e 



dwidqi / du)2dc\2 



1 



I9i| 



192 I 



7kil + kiP 7^2! + |g2|^ 



- t^i) - VT,{k^ + qix) - N{ky + qiy)^/2m 

1 



- UJ2) - Vp{k^ + q2x) - N{ky + q2yY/2m 



Y,{uJrn ~ uji~ UJ2) - Vp{kx + qix + q2x) - N{ky + qiy + q2yY I2m 
Integrating over qix and q2x and rescaling frequencies wi, L02 by Um and momenta qiy, q2y by {'^ijJi^2Y^^ ^ we obtain 

dqidq2\qiq2\ d^i dQ2 



Y^2{k,LUm) 



^\ 2/3 1/3 

2^ < -0 



(l + k-l|3)(l-|-|g-2|3)7„ f7f 
1 



(B2) 



2\/3A^gig2(f^if^2)i/=^ + 



+ r!i)2/3 + (1 _ r!2)2/3 + (f^^ + - - 1] 



where 57 1.; 



and (71,; 



g(1.2)!; 



r 



a. Large N expansion 

The two- loop self-energy to leading order in 1/iV is 
obtained by expanding the denominator in (|B2p to order 
1/N'^. The two momentum integrals are logarithmical 
with the lower limit set by 1/A^. We have 



T,2{k,LUrn) = - 



1 rff^i 



1/3 2/3 



(1 - r!i)2/3 + (1 - 172)'/' 



dn2 



-{^1+^2 - 1) 



2/3 



(B3) 



Integrating over fii and U2 in Eq. (|B3[) we obtain 
Eq. (|3.2p of the main text. 



6. Expansion near the mass-shell 

Eq. (|B3[) is valid when In TV is large and the result, 
Eq. (|3.2p . contains only the term linear in A. For small 



N, the momentum integration has to be done more ac- 
curately. Below we expand E2 near the mass shell and 
show that this expansion contains a non-analytic A^/^ 
term. This non-analytic term is comes from low-energy 
fermions and is enhanced by vertex corrections. 

We assume and then verify that the non-analytic term 
comes from the regions ili <C 1, 1 — 572 ^ 1 and fli <C 1, 
1 — f72 ^ 1. The contributions from these two regions 
are equal and we only focus on the contribution from 
rii <C 1, 1 - ^2 < 1- We introduce 17^ = 1 - O2 and 
make use of the identity 



dq2\q2\ 



1 



-00 1 + kiP 7-00 1 + k2P A + iBqiq2 
A 

54(^6 + B6) 



(B4) 



+ 108B In {B /A) \ ttA^ + ln(B /A) 



where in our case 

B = 2\/3iV(^^l)^/^ A = A- 



(f7i - Vt 



I \2/3 



(B5) 
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The source for non-analyticity is the term with In^ (B/A), 
and keeping only this term we obtain 



Introducing new integration variables, x = f^i/A'^/^, y 
nyA^/^, we re- write Eq. ([XTB)) as 



2ti 



ill 



dVL' 



, n2/3 



,1/3 



1/3 



X In^ 



1/3 



2/3 



(B6) 



E2(fc,u;„ 



■ 2/3 1/3 



2/3 



In 



1 



In 



,1/3 



1 + y2/3 + (2; - j/)2/3^ 



(B7) 



This expression contains contributions confined to the 
upper limit of the integration over x, but also contains 
universal contributions which come from x,y — 0{1) and 
is independent on the upper limit of integration. One can 
easily verify that the largest contribution of this kind 
comes from the cross-product of the two logarithmic fac- 
tors in the last bracket. Evaluating the integrals explic- 
itly we obtain Eq. (|3.16p of the main text. 



the two diagrams presented in Fig. [2l3. Compared to 
the two-loop diagram, these two diagrams describe cor- 
rections to the two spin-boson vertices which involve a 
boson with a frequency u)2 ~ w. Both diagrams give 
equal contributions, and we consider only the one from 
Fig. [Jb- The analytical expression for this diagram is 



2. Three-loop diagram 

We verified that the most singular contributions to 
fermionic self-energy at the three-loop level come from 



f duji dqi f duj2 dq2 f dwa dqs 



i,'^i)x(q2,t^2)x(q3,t^3)G(k + qi,Um - i^l) 



X G(k + qi + q2, ujm - t^i - W2)G(k + q2, - i^2)G(k + q2 + qs, - W2 - W3)G(k + qg, uj„i - ^408) 
duJidqi j duj2dq2 I dcusdqs- 



9 



(2^) 



193 I 



7|wi| + \qi\'^ 7|W2| + 1921'' 7l^^3| + ml 



1 



1 



T,{uj„i - t^i) - Vp{kx + qi^j;) - N{ky + qiyY /2m 
1 

t.{ujni - W3) - Vp{kx + qzx) - N{ky + qzyY/2m 

1 

i;(a;„i -ui- UJ2) - Vp{kx + qi^ + q2x) - N{ky + qiy + q2yY/2m 

1 

i;(w„i - UJ2- t^3) - v^{kx + q2x + qzx) - N{ky + q2y + q^yY /2m 



I](w,„ - U!2) - Vp{kx + q2x) ~ N(ky + q2y)'^/2r 



(B9) 
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Performing the integratfon over variables qix, q2x, and q^x we obtain 

3 



dqxdq2dq-i\qxq2qz\ 



_^ (1 + |g-i|3)(l + |g-2|3)(l + |g-3|3) ^1/3 J^_^^ ^1/3 J^_^^ ^ 
1 

2ViNq,q2in,n2y^^ + + ^ [(i - n,)^^ + (i _ ^12)2/3 + 1)2/3 _ 1] 

1 

2V3A^g293(f^2f^3)i/' +iA + i[{l- ris)^/^ + (1 - ^3)^/3 + (^ia + ^3 ~ l)^/^ - l] ' 



1/3 

3 



(BIO) 



where ili 2 



and gi 2 



9(1, 2)y 



Yjj. Compared to the two- loop diagram given by Eq. p.7p this one contains an 



additional integral over fl^ and additional denominator. 

We are interested in the corrections to the universal, non-analytic term in the self-energy. Accordingly, we introduce, 
as before, ri2 = 1 — and set Slj, Sli <C 1. We then obtain 



2ti 



dqidq2dq3\qiq2q3\ C ijC)' C '^^^ C ^^3 

(1 + k-i|3)(l + |g-2|3)(l + |g-3|3) ^ J^, ^1/3 J^, ^1/3 



(Bll) 



1 



-{2V5i)Nqiq2ny^ + A+ ^2^''^ + (f^i - ^)2/3 - {2V3 i) Nq2q3n'/^ + A + ^'^^^'^ + {^3 - il'^)^/^ 



y 2/3 



1/3 



y 2/3 



a. N = 

Integrating over dimensionless momenta qi and 52 in Eq. (jBlip and setting = we obtain 

I]3(fc,Wm) = ^ i Ujt{'^UJn'^ I / dQ 



™ ""0 / 1/3 , 
Jo 

2 dn^ 



1 



3 Ju'^ A + + (O3 - r!^)2/3 



(B12) 



The first line is the expression for one of the two two-loop 
self-energy diagram, the second line is the extra piece 
which represents the vertex correction. We see that the 
integral over is logarithmical. Using the fact that typ- 
ical ~ A3/2, we obtain, with logarithmical accuracy. 



S3(fc,tJm) = E2(fc,Wm)ln 



(B13) 



APPENDIX C: FERMIONS AT THE 
HALF-FILLED LANDAU LEVEL 



In this Appendix we give details of the calculations 
for a model of fermions interacting with a bosonic field 
whose static propagator scales as l/|g|. 



This is the result that we cited in the main text. 



1. Real-space propagator for N ^ 1 



b. Ftmte N 

At a finite N and A ^ N'^, the logarithmic divergence 
of the integral over fl^ in (|B12[) is cut by N'^ instead of 
A, and we have 



T,3{k,uj„i) = S2(fc,Wm)ln-^. (B14) 



The derivation of the form of G(r, t) is analogous to the 
derivation of Eq. (lASp . Using the self-energy = 
iXuJm In (co/'^m) and converting it to real frequencies, we 
obtain, with logarithmical accuracy. 



dk exp (ikr — iujt) 



(27r)2 \ujln{eo/uj) - tt 



(CI) 
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Introducing new variables x = Aw In (eo/i^) /cfe, y 
^kf/v-p and using the fact that 



, X h,, 6 cos a — sin a , „ 
cos a + y e-''^ = — , 6 > 0, C2 

1 + ^ 



we arrive at Eq. (|4.6p . 



2. Two-loop diagram 



The analytical expression for Y,2{<li^m) has the same 
structure as Eq. (|Bip and after the integration over 
x— component of momenta reduces to 



S2(fc,a;„j) = -A / duji I duo2 
Jo <y ujjj^—uji 



dqidq2\qiq2\ 



(7W1 + |gi|2)(7W2 + |g2p) 
1 



{Nqiq2/m) - efe + - wi) + - ^2) + + ^2 - w^) 

where 7 = j/{mu) — AirXm and gi_2 are y components of running momenta bounded by |gi,2| < fcp- 



(C3) 



a. Large N expansion 

As in the previous case, we expand the denominator in Eq. (jC3p in powers of [l/N) ^ 1 and cut the logarithms 
by 1/A^. To leading order in 1/A^, we have 



S2(fc,CJ„) 



Using 



and 



" duji f ™ duj2 I , ., 
/ < -efe + «A 

"^1 Jui^-uix ^2 I 



+ (cj„i - CJ2) In 



— W2 



+ (Wi + 0^2 - t^m) In 



Co 



UJi+ UJ2 — LOn 



" duJi 

^1 Juj^-o; 



C?ti;2 
1 0^2 6 



" dwi f'^™ duj2 

Jui^-uJi ^2 



-a;i)ln 

+ (wi + 0^2 - tt^m) In 



+ {uj,n - W2) In 



— w,„ln(eo/w„), 




^ + a;2 — ti-"' 

and substituting Eqs. (|C5p and (|C6p into Eq. (|3.2p we obtain Eq. (|4.4p of the main text. 



(C4) 



(C5) 



(C6) 



h. Universal term in the expansion near the mass-shell 

As in Appendix [B] we introduce new dimensionless variables qi = qi/ and 92 = 92 / \/j(^2- The integration 
over qi and 92 in Eq. (jCSp is now confined to |gi,2| < (eo/1^1,2)^^^, and Eq. (jCSp becomes, after proper rescaling 



V ^i, ^ ^ /"^^o /" ^^^il^il M92|g2l 

E2(K,t^m) = «A w,„ / dill / di22 / , -2 / , , -2 X 

1 



(C7) 



'iNa^/TlJhqiq2 + Aln (eo/w,„) + |(1 - 



In 



+ (l-r!2)ln 



l-f22 



(r!i + r!2-l)ln 



where i7i^2 = (^1,2/^™), a — {1/2'k){vf /u) and A is given by Eq. 
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As before, we expand E2 near the mass-shell in powers of A and search for a universal, non-analytical contributions 
from f2i ~ 0, 1 and vice versa. Restricting with the contribution from the first region and introducing 

= 1 — ^ 1, we obtain 

V ^i, \ x /"\n' f dqi\qi\ f dq2\q2\ 

2^2(k,i^m) = «Aw„ / dn^ / dill / / — — T 

7 1 + 91 7 1 + 92 

1 

^ 

-i7VaVnrgig2 + A In (eo/w™) + In + (f^i - In (ht^) 

(C8) 



The integrals over qi and 52 are confined to jqil < Qi = 
(eoM„)'/' (l/fli)i/2 and jgsl < Q2 = (eo/c^™)'^', re- 
spectively, the upper limits of the integrals over Qi and 
should not matter in this approximation. 
A simple experimentation shows that the universal, 
non-analytic contribution to E2 appears already at first 
order in A. Expanding (|C8p in A and integrating over 
qi and 52 we obtain with the logarithmic accuracy 

rl /.Oi 

In(lM) 

^ln(j.)+(f7i-fiyin(^)]' 

~ 2 1 
= -zAw„iAln (eo/wm) Inln — . (C9) 

"0 



The lower cutoff f2o is defined by rioln(l/rio) ^ 
Alneo/t^m at TV = and by f^o ln^(l/f^o) iV^ In'' TV 
at iV^ ^ A (up to extra logarithms). This leads to 
Eqs. (|iTT|) and fiTT)) in the main text. 
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